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The physical interpretation of time-resolved photoabsorption experiments is not as straightforward as for the 
more conventional photoabsorption experiments conducted on equilibrium systems. In fact, the relation between 
the transient photoabsorption spectrum and the properties of the examined sample can be rather intricate since 
the former is a complicated functional of both the driving pump and the feeble probe fields. In this work we 
critically review the derivation of the time-resolved photoabsorption spectrum in terms of the nonequilibrium 
dipole response function \ and assess its domain of validity. We then analyze \ in detail and discuss a few exact 
properties useful to interpret the transient spectrum during (overlapping regime) and after (nonoverlapping 
regime) the action of the pump. The nonoverlapping regime is the simplest to address. The absorption energies 
are indeed independent of the delay between the pump and probe pulses and hence the transient spectrum can 
change only by a rearrangement of the spectral weights. We give a close expression of these spectral weights 
in two limiting cases (ultrashort and everlasting monochromatic probes) and highlight their strong dependence 
on coherence and probe-envelope. In the overlapping regime we obtain a Lehmann-like representation of \ in 
terms of light-dressed states and provide a unifying framework of various well known effects in pump-driven 
systems. We also show the emergence of spectral sub-structures due to the finite duration of the pump pulse. 

PACS numbers: 78.47,jb,32.70.-n,42.50.Hz,42.50.Gy 


I. INTRODUCTION 

Filming a “movie” with electrons and nuclei as actors may 
sound fantasy-like, but it is de facto a common practice in 
physics and chemistry modern laboratories. With the im¬ 
pressive march of advances in laser technology, ultrashort 
(down to the sub-fs time-scale), intense (> 10 15 W/cm 2 ) 
and focussed pulses of designable shape, hereafter named 
pumps , are available to move electrons in real or energy 
space. By recording the photoemission or photoabsorption 
spectrum produced by a second, weak pulse, hereafter named 
probe, impacting the sample with a tunable delay from the 
pump a large variety of ultrafast physical and chemical pro¬ 
cesses can be documented. 1-6 Time-resolved pump-and-probe 
(P&P) spectroscopies have revealed the formation and dy¬ 
namics of excitons, 7-11 charge-transfer excitations, 12-20 auto- 
ionized states 21,22 and light-dressed states, 23-37 the evolution 
of Fano resonances, 38-43 the screening build-up of charged 
excitations, 44 ^ 47 the transient transparency of solids, 4849 the 
motion of valence electrons, 50-55 the band-gap renormaliza¬ 
tion of excited semiconductors, 51 ^ 58 how chemical bonds 
break 59-62 and other fundamental phenomena. 

A suited P&P spectroscopy to investigate charge-neutral 
excitations is the time-resolved (TR) photoabsorbtion (PA) 
spectroscopy. 63-65 It is well established that PA spectra of 
equilibrium systems are proportional to the dipole-dipole re¬ 
sponse function x, 66-69 an extremely useful quantity to under¬ 
stand and interpret the experimental results. In pump-driven 
systems the derivation of a mathematical quantity to interpret 
TR-PA spectra is slightly more delicate and, in fact, several 
recent works have been devoted to this subject. 7(1-74 The diffi¬ 
culty in constructing a solid and general TR-PA spectroscopy 
framework (valid for general P&P envelops, durations and de¬ 
lays and for samples of any thickness) stems from the fact that 


the probed systems evolve in a strong time-dependent electro¬ 
magnetic (em) field and hence (i) low-order perturbation the¬ 
ory in the pump intensity may not be sufficiently accurate and 
(ii) separating the total energy per unit frequency absorbed by 
the system into a pump and probe contribution is questionable. 
Furthermore, due to the lack of time-translational invariance 
the TR-PA spectrum is not an intrinsic property of the pump- 
driven system, depending it on the shape of the probe field 
too. 

We can distinguish two different approaches to derive a 
TR-PA formula: the energy approach, 7() ' 72 ' 74 which aims at 
calculating the energy absorbed from only the probe, and the 
Maxwell approach, 71,73 ' 75 which aims at calculating the trans¬ 
mitted probe field (these approaches are equivalent for opti¬ 
cally thin and equilibrium samples). We carefully revisit the 
energy approach, highlight its limitations and infer that it is 
not suited to perform a spectral decomposition of the absorbed 
energy. We also re-examine the Maxwell approach and pro¬ 
vide a derivation of the TR-PA spectrum in non-magnetic sys¬ 
tems without the need of frequently made assumptions like, 
e.g., slowly-varying probe envelops or ultrathin samples. The 
final result is that the TR-PA spectrum can be calculated from 
the single and double convolution of the nonequilibrium re¬ 
sponse function \ with the probe field. 

For the physical interpretation of TR-PA spectral features a 
Lehmann-like representation of the nonequilibrium \ would 
be highly valuable, as it is in PA spectroscopy of equilibrium 
systems. In this work we discuss some exact properties of 
the nonequilibrium \ and of its convolution with the probe 
field. When the probe acts after the pump (nonoverlapping 
regime) \ can be written as the average over a nonstationary 
state of the dipole operator-correlator evolving with the equi¬ 
librium Hamiltonian of the sample. In this regime the TR- 
PA spectrum is nonvanishing when the frequency matches the 
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difference of two excited-state energies. As these energies are 
independent of the delay r between the pump and probe field, 
only the spectral weights can change with r (not the absorp¬ 
tion regions). 76 We discuss in detail how the spectral weights 
are affected by the coherence between nondegenerate excita¬ 
tions and by the shape of the probe field. A close expression 
is given in the two limiting cases of ultrashort and everlasting 
monochromatic probes. 

The overlapping regime is, in general, much more compli¬ 
cated to address. The absorption energies cease to be an in¬ 
trinsic property of the unperturbed system and acquire a de¬ 
pendence on the delay. Nevertheless, an analytic treatment 
is still possible in some relevant situations. For many-cycle 
pump fields of duration longer than the typical dipole relax¬ 
ation time, we show that a Lehmann-like representation of 
the nonequilibrium \ in terms of light-dressed states can be 
used to interpret the TR-PA spectrum. We provide a uni¬ 
fying framework of well known effects in pump-driven sys¬ 
tems like, e.g., the AC Stark shift, the Autler-Townes splitting 
and the Mollow triplet. More analytic results can be found 
for samples described by a few level systems. In this case, 
from the exact solution of the nonequilibrium response func¬ 
tion with pump fields of finite duration we obtain the dipole 
moment induced by ultrashort probe fields. The analytic ex¬ 
pression shows that (i) the r-dependent renormalization of the 
absorption energies follows closely the pump envelope and (ii) 
a spectral sub-structure characterized by extra absorption en¬ 
ergies emerges. 

The paper is organized as follows. In Section II we briefly 
review the principles of TR-PA spectroscopy measurements. 
We critically discuss the energy approach to PA spectroscopy 
in equilibrium systems and highlight the limitations which 
hinder a generalization to nonequilibrium situations. The con¬ 
ceptual problems of the energy approach are overcome by 
the Maxwell approach which is re-examined and used to de¬ 
rive the transmitted probe field emerging from non-magnetic 
samples of arbitrary thickness, without any assumption on the 
shape of the incident pulse. The nonequilibrium dipole-dipole 
response function \ is introduced in Section III and related 
to the transmitted probe field. We analyze x i n the nonover¬ 
lapping regime in Section IV and in the overlapping regime 
in Section V. Finally, a class of exact solutions in few-level 
systems for overlapping pump and probe fields is presented 
in Section VI. Summary and conclusions are drawn in Sec¬ 
tion VII. 


H. TIME-RESOLVED PHOTOABSORPTION 
SPECTROSCOPY 


In this section we briefly revisit the principles of PA for sys¬ 
tems in equilibrium, and subsequently generalize the discus¬ 
sion to the more recent TR-PA for systems driven away from 
equilibrium. The aim of this preliminary section is to high¬ 
light the underlying assumptions of orthodox equilibrium PA 
theories, identify the new physical ingredients that a nonequi¬ 
librium PA theory should incorporate, and eventually obtain 
a formula for the TR spectrum which is a functional of both 
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FIG. 1: (Color online) Illustration of a PA experiment. 


the pump and probe fields. Except for a critical review of the 
literature no original results are present in this section. 


A. Experimental measurement 


Consider a system in equilibrium and irradiate it with some 
feeble light (perturbative probe). In Fig. 1 we show a snapshot 
at time t of a typical equilibrium PA experiment. The incident 
light is described by the electric and magnetic fields e(t) and 
b(t) (left side of the sample) whereas the transmitted light is 
described by the em fields e'(f), b'(f) (right side of the sam¬ 
ple). The experiment measures the total transmitted energy 
E'. This quantity is given by the energy flow (or equivalently 
the Poynting vector) integrated over time (the duration of the 
experiment) and surface. Denoting by S the cross section of 
the incident beam we have (here and in the following integrals 
with no upper and lower limits go from —oo to +oo) 

E' = S-^~ [dt |e'(i) x b'(f)|. (1) 

47T J 

The integral in Eq. (1) is finite since the em fields used in an 
experiment vanish outside a certain time interval. In vacuum 
the electric and magnetic fields are perpendicular to each other 
and their cross product is parallel to the direction of propa¬ 
gation. Taking into account that |b'| = |e'| the transmitted 
energy in Eq. (1) simplifies to 

E'= S^ J dt\e'{t)\ 2 . (2) 


From Eq. (2) we see that the transmitted energy E' depends 
on the temporal shape of the electric field. This dependence 
can be exploited to extract the energy of the neutral excitations 
of the sample. A typical, systematic way of varying the tem¬ 
poral shape consists in probing the sample with monochro¬ 
matic light of varying frequency. Taking into account that the 
electric field is real, its Fourier transform reads 


r OO 7 

= l ^-e’(ut)e- iut + c.c. (3) 

where c.c stands for “complex conjugate”. Unless otherwise 
defined, quantities with the tilde symbol on top denote the 
Fourier transform of the corresponding time-dependent quan¬ 
tities. Inserting Eq. (3) back into Eq. (2) we find 


fSleW- 
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For monochromatic light of frequency u>q (in the time interval 
of the experiment) the transmitted field e' (ut) is peaked at w ~ 
wo and hence E' ~ |e'(co»o)|where Aw is the width 

of the peaked function e'(w). Therefore the quantity 

T^ , (w)=5^|g'(w)| 2 , w>0 (5) 

can be interpreted as the transmitted energy per unit fre¬ 
quency. 

Alternatively W'{ w) could be measured using fields of ar¬ 
bitrary temporal shape and a spectrometer. As this is also the 
technique in TR-PA experiments, and the method to elaborate 
the data of the real-time simulations of Section VI is based on 
this technique, we shortly describe its principles. The spec¬ 
trometer, placed between the sample and the detector, splits 
the transmitted beam into two halves and generates a tunable 
delay 5 for one of the halves. The resulting electric field at the 
detector is therefore i (e'(f) + e'(t — <$)), and the measured 
transmitted energy is 


E'(S) =S±Jdt 


[ dt 

e'(f) + e'(t — 8) 

to 

1 


( 6 ) 


The PA experiment is repeated for different delays <5, and the 
results are collected to perform a cosine transform 


E'{y) 


J dS E'{8) cos (v8). 


(7) 


The relation between E' and W is readily found. Using Eq. 
(3) we get 




2 1 + COS (W(5) 


( 8 ) 


Inserting this result into Eq. (7) and taking into account the 
identity J dS cos(w5) cos(iz£) = tt [<5(w + v) + <5(w — v)\ we 
find E'{v) = 7r 8{v)E' + <S^r |e'(i/)| 2 . Thus for every v ^ 0 
we have W'{v) = AE'(v). 

The PA experiment can be repeated without the sample to 
measure the energy per unit frequency W(ut) of the incident 
beam. The difference 


5(w) = W(w) - W'(w) > 0 (9) 


is therefore the missing energy per unit frequency. How to 
relate this experimental quantity to the excited energies and 
excited states of the sample is well established and will be 
reviewed in the next Section. 

The main novelty introduced by TR-PA experiments con¬ 
sists in probing the sample in a nonstationary (possibly driven) 
state. The sample is driven out of equilibrium by an intense 
laser pulse described by the em fields E(f) and B(t), and sub¬ 
sequently probed with the em fields e(t) and b(t), see Fig. 2. 
We refer to E (t) and B(f) as the pump fields. To extract in¬ 
formation on the missing probe energy the transmitted pump 
field is not measured, see again Fig. 2. As the sample is not in 
its ground state the transmitted beam is also made of photons 
produced by the stimulated emission. These photons have the 
same frequency and direction of the probe photons and, there¬ 
fore, the inequality in Eq. (9) is no longer guaranteed. 



FIG. 2: (Color online) Illustration of a TR-PA experiment. 


B. Energy approach: PA in equilibrium 

In this section we obtain an expression for the spectrum 
S(ui) of systems initially in equilibrium in their ground state 
(the finite-temperature generalization is straightforward). We 
use an approach based on the energy dissipated by the sam¬ 
ple and highlight those parts in the derivation where the hy¬ 
pothesis of initial equilibrium and weak em field is used. We 
advance that this approach cannot be generalized to nonequi¬ 
librium situations. 

For a system driven out of equilibrium by an external (trans¬ 
verse) electric field E ext the energy absorbed per unit time, 
i.e., the power dissipated by the system, is 

V(t) = -jdr J(rf) • E ext (rf), (10) 

where J is the current density and e = —1 is the electric 
charge. This well-known formula is valid only provided that 
the electric field generated by the induced current J is much 
smaller than E ext . For the time being let us assume that 
this is the case. We also assume that the probed systems are 
nanoscale samples like atoms and molecules or thin slabs of 
solids. Then the wavelength of the incident em field is typi¬ 
cally much larger than the longitudinal dimension of the sam¬ 
ple and the spatial dependence of E ext can be ignored. Writ¬ 
ing J = V(r • J) — (V • J)r, discarding the total divergence 
and using the continuity equation the power becomes 

Vit) = - J dr r • E ext (f), (11) 

where n is the electron density. The integral of the power be¬ 
tween any two times t\ and t 2 yields the difference between 
the energy of the system at time t 2 and the energy of the sys¬ 
tem at time tp 


E sys (t 2 ) - E sys (ti) = - dt E ext (f) • (12) 

where we found convenient to define the dipole moment 

d(f) = J drrn(ri). (13) 

We observe that nowhere the assumption of small external 
em fields and/or the assumption of a system in equilibrium 
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are made in the derivation of Eq. (12). It is also important 
to emphasize the semiclassical nature of Eq. (12). Suppose 
that the sample is initially in its ground state with energy E g . 
We switch the em field on at a time t = t on and switch it 
off at a time t = f Q g. According to Eq. (12) the difference 
E aya (t) — Eg remains constant for any t > f 0 ff- Physically, 
however, this is not what happens. At times t S> f 0 ff the sam¬ 
ple is back in its ground state since there has been enough 
time to relax (via the spontaneous emission of light). Hence 
the correct physical result should be E sys (t —>• oo) — E g = 0. 
In Eq. (12) the description of the em field is purely classical 
and does not capture the phenomenon of spontaneous emis¬ 
sion. Nevertheless, a spontaneous emission process occurs on 
a time scale much longer than the duration of a typical PA ex¬ 
periment. The semiclassical formula is therefore accurate for 
times t > t 0 e and, consequently, the quantity 

f d 

E ahs = ~ dt E ext (f) • — d(t) (14) 

can be identified with the increase in the energy of the sample 
just after the em field has been switched off. We refer to E abs 
as the absorbed energy. Let us see how to relate Eq. (14) to 
the missing energy S(u>) measured in an experiment. 

In equilibrium PA experiments E oxt = e is the probing 
field discussed in Section II A. Let E = f (w) and E' = 
f be the total energy of the incident and transmitted 

beam respectively. Then the difference E — E’ is the energy 
transferred to the sample, i.e., the absorbed energy of Eq. (14) 

E-E , = £ abs . (15) 

We write the dipole moment d = d cq + d ; , as the the sum of 
the equilibrium value d eq and the probe induced variation d p . 
Since d eq is constant in time the absorbed energy in frequency 
space reads 

/•OO 

-ftabs = i / — W e*(cj) • d p (v) + c.c. (16) 

Jo ^ 

Taking into account Eq. (15) and the definition of S in Eq. (9) 
we also have 

f d 

Eabs = / —S(cj). (17) 

We now show that the r.h.s. ofEqs. (16) and (17) are the same 
because the integrands are the same. The transmitted em field 
e'(uj) at frequency uj depends, to lowest order in e, only on 
e(w) at the same frequency uj since the system is initially in 
equilibrium (hence invariant under time translations). This 
implies that if the probe field has N frequencies uj±, ... ,u>n 
then the total missing energy E — E l is the sum of the missing 
energies of N independent PA experiments carried out with 
monochromatic beams of frequencies u> i, • • •, ujn- The same 
is true for the energy absorbed by the sample: d p (w) depends 
only on e(w) since the probe-induced dipole moment is linear 
in e. Therefore for systems in equilibrium and to lowest order 
in the probing fields we can write 

(18) 


The approaches to calculate the right hand side of Eq. (18) 
can be grouped into two classes. In one (recently emerg¬ 
ing) class one perturbs the system with an em field e(t), 
calculates the time-dependent dipole moment either by solv¬ 
ing the Schrodinger/Liouville equation or by using other 
methods, 77-87 and then Fourier transforms it. The other (more 
traditional) class avoids time-propagations and works directly 
in frequency space. To lowest order in e the Kubo formula 
gives 

d pA t ) = '^2 [ (19) 

3 

where x is the (retarded) dipole-dipole response function. For 
a system with Hamiltonian H in the ground state 'f (/ ) of en¬ 
ergy Eg we have 

= 0(t-t')(^ g \e l ^ t d i e~ l ^ < ' t ~ t ' , d j e~ l ^ t -H.c-IT'p) 
= 6{t - t')^ g \d i e~ i ^~ Ea ' , ^ t ~ t '' > dj - H.c.|$ 9 ) (20) 

where H.c. stands for “hermitian conjugate” and d, is the i- 
th component of the dipole-moment operator. As expected 
the equilibrium response function x depends on the time- 
difference only. Fourier transforming Eq. (19) and inserting 
the result into Eq. (18) we get 

S{u) =a;y^g*(a;)AjMej(a;), (21) 

ij 

with Cij{uj) = i [xij(oj) — Xjj(w)]. The response function 
x(w) can be calculated by several means without performing 
a time propagation. From the Lehmann representation of x it 
is easy to verify that C is positive semidefinite for positive fre¬ 
quencies and negative semidefinite otherwise. Consequently 
S is manifestly positive, in agreement with Eq. (9). 

C. Energy approach: PA out of equilibrium 

In a typical TR-PA experiment both the pump and probe 
fields are very short (fs-as) laser pulses with a delay r be¬ 
tween them. If r < 0 then the probe acts before the pump 
and we recover the PA spectra of equilibrium systems. On the 
other hand if r > 0 then S(ui) acquires a dependence on r. 

This dependence can be used to follow the evolution of the 
system in real time. However, for the physical interpretation 
of what we are actually following it is necessary to generalize 
the equilibrium PA theory to nonequilibrium situations. 

Let the external electric field E ext be the sum of the pump 
field E and probe field e, i.e., E ext = E + e. In this case Eq. 

(14) yields the total energy absorbed by the system. As the 
experiment detects only the energy of the transmitted probe 
field the use of the energy approach for TR-PA is not straight¬ 
forward. One might argue that the energy absorbed from the 
probe is given by Eq. (14) in which E ex t —>• e 

/ d 
dt e(t) ■ ^-< 1 ( 7 ). 


S(ui) = iui e*(w) ■ d p (u;) + c.c. 


( 22 ) 
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However, this formula cannot be always correct. Suppose that 
the pump field is also feeble and can be treated as a small per¬ 
turbation. Then, the transmitted probe field e'(w) depends on 
E(w) + e(w) (linear response theory). These fields are in¬ 
dependent of position inside the sample. In a larger space, 
like that of the laboratory, they do depend on r and this de¬ 
pendence specifies the direction of propagation. Let E(Ku) 
and e(kw) be the spatial Fourier transform of the pump and 
probe fields. For isotropic systems e'(kcj) depends only on 
e(kcu) since E(kcu) vanishes for k parallel to the direction of 
propagation of the probe. This implies that the missing energy 
per unit frequency is independent of the pump, a conclusion 
which is not in agreement with Eq. (22). In fact, E a b s de¬ 
pends on E whenever the pump-induced variation of d is not 
orthogonal to e. 

To cure this problem we could write d = dp + d p , where 
dp is the value of the dipole moment when only the pump 
field is present whereas d p is the probe-induced variation, and 
say that the missing probe energy is 

E abs = - J dt e(t) ■ ^d p (f). (23) 

This expression is by construction correct for perturbative 
pumps. For pumps of arbitrary strength Eq. (23) cannot be 
proved or disproved using exclusively energy considerations. 
For the sake of the argument, however, let us assume that Eq. 
(23) is the correct missing energy. There is still a conceptual 
problem to overcome if we are interested in the missing en¬ 
ergy per unit frequency. For strong pump fields the sample is 
in a nonstationary state and hence the transmitted probe field 
e'(oj) depends on the entire function e, not only on the value 
e(w) at the same frequency. Thus the reasoning made below 
Eq. (17) does not apply. In particular if e is monochromatic 
then e' (as well as d p ) is, in general, not monochromatic. Con¬ 
sequently S oc |e| 2 — |e'| 2 is, in general, not monochromatic 
either. If we used the formula in Eq. (18) we would instead 
find that S is peaked at only one frequency since e is peaked 
at only one frequency. To overcome these problems one has 
to abandon the energy approach and calculate explicitly the 
transmitted probe field. 

D. Maxwell approach 


unit vectors along the propagation direction of the pump and 
probe fields respectively. In TR-PA experiments these vectors 
are not parallel for otherwise the detector would measure the 
transmitted pump intensity too. The time-dependence of the 
macroscopic current density arises when the pump and probe 
fields interact with the electrons in the sample. For trans¬ 
verse pump and probe fields and for isotropic systems (J) is 
the sum of transverse waves propagating along the directions 
QKN + qkn with Q and q integers, and K and k the pump 
and probe wave numbers respectively. Consequently, the to¬ 
tal electric field too is the sum of waves propagating along 
QK'N + qkh, and Eq. (24) can be solved for each direction 
separately. 

We define £ p and (J) p as the wave of the electric field and 
current density propagating toward the detector (hence Q = 0 
and q = 1). The vectors £ p and (J) p depend on the spatial po¬ 
sition r only through x = n ■ r (transverse fields) and are par¬ 
allel to some unit vector e p lying on the plane orthogonal to n 
(the generalization to multiple polarization is straightforward, 
see Section VIB): £ p = e p £ p (xt) and (J) p = e p J p (xt). 
Equation (24) implies that 

d 2 £ p 1 d 2 £ p 4-7T dJ p 

8x 2 c 2 dt 2 c 2 dt’ C ^ 

which is a one-dimensional wave equation that can be solved 
exactly without assuming slowly-varying probe envelops 71 
or ultrathin samples. 73 The electric field £ p is the sum of 
an arbitrary solution h(t — x/c) of the homogeneous equa¬ 
tion and of an arbitrary special solution s(xt): £ p (xt) = 
h(t — x/c) + s(xt). Without loss of generality we take the 
boundaries of the sample at x = 0 and x = L. Let e(t — x/c ) 
be the amplitude of the incident probe field which at time 
t = 0 is localized somewhere on the left of the sample. Im¬ 
posing the boundary condition £ p (x0) = e{—x/c) we then 
obtain 

£ p (xt) = e(t — x/c) — s(x — ct 0) + s(xt). (26) 

The special solution s(xt) is found by inverting the one¬ 
dimensional d’Alambertian □ = — -pz-jyp- The Green’s 

function G(xt) solution of DG(x t) = S(x)6(t) is 

G(xt) = ~^0(t)x[-ct,ct](x), (27) 


To overcome the difficulties of the energy approach we use 
the Maxwell equations to calculate explicitly the transmitted 
probe field. In a nonmagnetic medium the total electric field 
£, i.e., the sum of the external and induced field, satisfies the 
equation 75 - 88 


where X[a,b] (x) = 1 if x £ (a, b) and zero otherwise. There¬ 
fore the special solution reads 


s(xt) 


2n 

c 





dx' 


' —x—cit—t') 


dJ p {x't') 

dt' 


(28) 


1 d 2 £ _ 4tt d(3) 
c 2 dt 2 c 2 dt 


(24) 


where (J) (rf) is the macroscopic current density, i.e., the spa¬ 
tial average of the current density over small volumes around 
r. In the derivation of Eq. (24) one uses that V x V x £ = 
V(V • £) — V 2 £ and that V • £ = 0 since the sample is 
charge neutral (in a macroscopic sense). Let N and n be the 


Without any loss of generality we can choose the time t = 0 
as the time before which nor the pump and neither the probe 
have reached the sample. Then J p = 0 for t < 0 and hence 
s(x 0) = 0 for all x. We conclude that 

£ p (x t) = e(t — x/c) + s(x t), 

with s given in Eq. (28). 


(29) 
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We are interested in the electric field on the right of the 
sample, i.e., in x = L, since this is the detected field. Let us 
therefore evaluate Eq. (28) in x — L. Taking into account that 
Jp(x't') is nonvanishing only for x' £ (0 ,L) and t' > 0 we 
have 


s{Lt) 


27r 
c 




dJ p (x't') 

dt' 



dr (J)p(rf), 


(30) 


m. NONEQUILIBRIUM RESPONSE FUNCTION 

From the definition in Eq. (9) the spectrum of a (equilibrium 
or nonequilibrium) PA experiment is given by 

S(w) =S ^ 0 g M| 2 “ le'MI 2 ) • (33) 

Since e = e p e and rl ? , = e p d p are both parallel to e p , so it 
is e' = e p e'. Then, the Fourier transform of Eq. (32) yields 
e'(w) = e(uj) d p (oj), and the spectrum in Eq. (33) can 

be rewritten as 


where in the second line we integrated over the volume 
V = SL of the sample. Using again the identity (J) p = 
V(r • (J) p ) — r(V • (J) p ) and extending the integral over all 
space (outside V the current density vanishes) we can rewrite 
Eq. (30) as s(Lt) = ~^ £ p • f dr r(V ■ (J) p ). Substituting 
this result into Eq (29) and taking into account the continuity 
equation V • (J) p = —d{n) p /dt, where (n) p is the macro¬ 
scopic probe-induced change of the electronic density propa¬ 
gating along n, we eventually obtain the transmitted electric 
field" 


£ p (Lt) = e(t) + |^£p • ^ J drr(n) p {rt), (31) 

where we discarded the delay L/c in the first term on the right 
hand side. The volume integral is the probe-induced dipole 
moment propagating along n. In general this is not the same 
as the full probe-induced dipole moment dp defined above 
Eq. (23) since, to lowest order in the probe field, d p is the 
sum of waves propagating along fcn + Q/\'N. Although it is 
reasonable to expect that the wave propagating along n (i.e., 
with Q = 0) has the largest amplitude, it is important to bear 
in mind this conceptual difference. In fact, in equilibrium PA 
experiments d p and f dr r (n) p are the same due to the ab¬ 
sence of the pump. For not introducing too many symbols we 
redefine d p = f dr r (n) p . Then, by definition, d p is parallel 
to s p and we can cast Eq. (31) in vector notation as 

e'(f) =e(f) + ^^dp(f). (32) 

Equation (32) relates the transmitted probe field to the 
quantum-mechanical average of the probe-induced dipole mo¬ 
ment, and it represents the fundamental bridge between theory 
and experiment. The result has been derived without assum¬ 
ing that the wavelength of the incident field is much larger 
than the longitudinal dimension of the sample (for thick sam¬ 
ples £p can be substantially different from e and the quantum 
electron dynamics should be coupled to the Maxwell equa¬ 
tions). Equation (32) can be used to calculate the missing 
energy per unit frequency of pump-driven systems. Notewor¬ 
thy Eq. (32) is valid for positive and negative delays between 
pump and probe, as well as for situations in which pump and 
probe overlap in time or even for more exotic situations in 
which the pump is entirely contained in the time-window of 
the probe. 


S(w) = — 21m e*(oj)dp(uj)^ — ^ wd p (w) 


(34) 


At the end of Section IIC we criticized the energy ap¬ 
proach since it predicts a single-peak spectrum for monochro¬ 
matic probes. Let us analyze Eq. (34) for the same case. 
For monochromatic probes of frequency oj q the first term in 
Eq. (34) vanishes for ui ojo whereas the second term 
is nonvanishing at the same frequencies of the transmitted 
probe field, see Eq. (32), in agreement with the discussion 
at the end of Section IIC. The quadratic term in the dipole 
moment is usually discarded in equilibrium PA calculations 
since d p and e oscillate at the same frequencies, and typically 
|e(w)| 3> / Sc)\u>d p {ui)\. If we discard the last term in Eq. 
(34) then we recover the spectrum of Eq. (18) of the energy 
approach. 

For the physical interpretation of nonequilibrium PA spec¬ 
tra it is crucial to understand the physics contained in the 
nonequilibrium dipole-dipole response function. In fact, d p 
can be calculated from the scalar version of the Kubo formula 
inEq. (19), i.e., 89 



dt'x{tfi')£ p (t'). 


(35) 


Here the scalar dipole-dipole response function \ is defined 
according to x = £ p,iXij £ p,j and £p * s tota l electric 
field (£ p ~ e if the induced field is small). Unfortunately, 
for pump-driven systems a Lehmann-like formula for x does 
not exist due to the presence of a strong time-dependent per¬ 
turbation in the Hamiltonian. By introducing the evolution 
operator U (t, if) from t' to t of the system without the probe, 
the nonequilibrium dipole-dipole response function reads 


QC(M') = 9(t-t') 

x g \U{t 0 ,t)dU(t,t')dU(t',t 0 ) -H.c.l’Tg) 

(36) 


where d = e p • d and / 0 is any time earlier than the switch-on 
time of the pump and probe fields. As anticipated the nonequi¬ 
librium x depends on t and t! separately. It is clear from Eq. 
(36) that x does not have a simple representation in terms of 
the many-body eigenstates and eigenenergies of the unper¬ 
turbed system. It is also easy to verify that Eq. (36) agrees 
with Eq. (20) in the absence of the pump. 

As a final remark before presenting some exact properties 
of x, we observe that in equilibrium, see Eq. (21), the ra¬ 
tio S{u})/\e{u>)\ 2 is independent of the probe field, i.e., it is an 
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intrinsic property of the sample. This is not true in nonequilib¬ 
rium, even if we discard the last term in Eq. (34) and approxi¬ 
mate £ ~ e. The physical interpretation of nonequilibrium PA 
spectra cannot leave out of consideration the shape and dura¬ 
tion of the probe, and the relative delay between pump and 
probe. In the next sections we discuss two relevant situations 
for interpreting the outcome of a TR-PA experiment. 


IV. NONOVERLAPPING PUMP AND PROBE 

Let us consider the case of a probe pulse acting after the 
pump pulse. We take the time origin t = 0 as the switch-on 
time of the probe. Then the pump acts at some time t = — r < 
0. For t > 0 the probe-induced variation of the dipole moment 
can be calculated from Eq. (35) with lower integration limit 
t' = 0. As we only need X for t, t' >0 we have U(t 0l t) = 
U(to,0)e l ^ t and similarly U(t', to) = e~ lEt U(0, to), with 
H the unpertubed Hamiltonian of the system. Defining T) = 
U(f) 1 fo)!^) as the quantum state of the system at time t = 0, 
the response function in Eq. (36) becomes 

i X (t,t') = 6(t - t')('H\e i6t de- iA{ - t - t ' ) de~ iflt ' -H.c.|’T) 

(37) 

which closely resembles the equilibrium response function of 
Eq. (20). In fact, without pump fields |'T) = e lEgt ° and 
Eq. (37) reduces to the equilibrium response function. We 
emphasize that in the presence of pump fields Eq. (37) is valid 
only for t, t' > 0. 

We expand the quantum state |'T) = E a c a |4r Q ) in terms 
of the many-body eigenstates |’T Q ,) of H with eigenenergy 
E a . The coefficients c a = c a e~ lEo,T depend on the delay 
r between the pump and the probe, c a being the expansion 
coefficients of the state of the system at the end of the pump. 
Inserting the expansion in Eq. (37) and using the complete¬ 
ness relation E 7 |'T 7 )('T 7 | = 1 we find 

ix(t, t') = 0(t - t') ^2 <£c /J e’ n “'» t e’ a '* t d a7 d 7 p - c.c. 

(38) 

where we defined the dipole matrix elements d aj = 
(SRq.|cZ|'I 7 7 ) and the energy differences fl a7 = E a — E 1 . In 
the following we use this result to study the outcome of a PA 
experiment in two limiting cases, i.e., an ultrashort probe and 
a monochromatic probe. 


A. Ultra-short probe 

The probe fields used in TR-PA experiments are ultrashort 
laser pulses. For optically thin samples the em field generated 
by the probe-induced dipole moment is negligible and does 
not substantially affect the quantum evolution of the system. 
Therefore, we can calculate d p from Eq. (35) with £ = e. For 
a delta-like probe e(t) = eo(5(f), hence e(w) = eo, we find 

d p [t) = —ie o ^2 c* C 0 e’ n °" 1 ,t d a 7 d 7( 8 + c.c. (39) 

a/3 7 


where we used the explicit form of the response function in 
Eq. (38). Fourier transforming this result we obtain the spec¬ 
trum 


sv) = 


X 


— 2 weg ^2 I m [e zn “ /5T c* a cpd OL1 d 1 ^ 
a/3 7 


u> — £i 7a + ir\ ui + fl 7j g + irj 


(40) 


where r) is a positive infinitesimal and we discarded the 
quadratic term in d p (thin samples). In Eq. (40) the depen¬ 
dence on the delay r enters exclusively through the phase fac¬ 
tors and, consequently, it is only responsible for modulating 
the amplitude of the absorption peaks. The position of the 
peaks is instead an intrinsic property of the unperturbed sys¬ 
tem. Thus, a change in the peak-position (discrete spectrum) 
or in the onset of a continuum (continuum spectrum) due to r 
should not be attributed to a change of the many-body ener¬ 
gies but to a redistribution of the spectral weights. 

Let us discuss Eq. (40) in some detail. For a system in 
equilibrium in the ground state (no pump) c Q = 1 for a = g 
and c a = 0 otherwise, and the spectrum reduces to 


S(u) = 27rweQ ^2 1 2 (<^( w ~ ^ 79 ) — d(cu + fl 7g )) ■ 

7 

(41) 

Since E 7 > E g the spectrum is nonnegative, in agreement 
with Eq. (9). In particular the height of the peak at some fre¬ 
quency w 0 is given by h g = 27r|w 0 |eg E 7: n 7? =| Wo J \ d 9 "t\ 2 ^ 
0. It is also interesting to consider the hypothetical situation of 
a pump pulse which brings the system from the ground state 
to an excited state | d/) = |4> x ) with energy E x . As the system 
is stationary this is the simplest example of a nonequilibrium 
situation. In the stationary case we have c a = 1 for a = x 
and c a = 0 otherwise, and hence the spectrum is again given 
by Eq. (41) with the only difference that the subscript “g" is 
replaced by the subscript “x”. Since E x is not the lowest en¬ 
ergy the positivity of the spectrum is no longer guaranteed. In 
fact, the height of the peak at frequency wq is 


h x = 27rw 0 eo ( ^2 I dx tI 2 ~ ^2 I d ^f) ( 42 ) 

-y:0 7 ;£=ctjQ wq 


which can be either positive or negative. The sign is positive if 
the absorption rate is larger than the rate for stimulated emis¬ 
sion and negative otherwise. We observe that in the stationary 
case the spectrum is independent of the delay. 

The most general situation is a system in a nonstationary 
state. From Eq. (40) the peak intensity at some frequency oj 0 
reads 


h = 27rwoeo EE± E Re [ e in °“> T c* cpd a7 d 7 p\ , 

a/3 ± 7:S7 7Q! =±aJo 

(43) 


where we introduced the short-hand notation E± 4 ± = 
A + +A_, with A an arbitrary mathematical expression. If T) 
is a superposition of degenerate eigenstates, hence c a ^ 0 for 
E a = E and c a = 0 otherwise, then the system is stationary 
and the height is independent of the delay. The dependence on 
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r is manifest only for |^) a superposition of nondegenerate 
eigenstates. The simplest example is a system in a superpo¬ 
sition of two eigenstates with energy E a and Eb, real coeffi¬ 
cients c a and ci, and real dipole matrix elements. In this case 
Eq. (43) yields h = h a + hb + h a b cos(fl a {,r) where h x=a ^ is 
defined as in Eq. (42) and 

hab — c 2tTUJq€.qCq i Ci ) EE ± E (I,, ', . 

a=a,b ± 7:Q 7a ==tu;o 

The spectral fingerprint of a nonstationary system is the mod¬ 
ulation of the peak intensities with r. These coherent oscilla¬ 
tions have been first observed in Ref. 54. 


B. Monochromatic probe 

The induced electric field of thick samples is not negligible 
and can last much longer than the external probe pulse. In 
this case the quantum evolution of the system should be cou¬ 
pled to the Maxwell equations to determine the total field self- 
consistently. 23 ' 35,40,70 ' 75 ’ 90 97 For a typical sub-as pulse e(t) 
centered around a resonant frequency wo the total electric field 
£(t) is dominated by oscillations of frequency wo decaying 
over the same time-scale of the induced dipole moment (in 
atomic gases this time-scale can be as long as hundreds of fs). 
Let us explore the outcome of a TR-PA experiment for a to¬ 
tal field of the form, e.g., £(t) = £oQ(t)sm(u>ot), wo > 0. 
Taking into account Eq. (38) we find 

gi(r2 aj a=bo;o)t _ 

-b C.C. 

±(fl 7j a ± wo) 

(44) 

If I'T) = |4/a;) (stationary system) then c a = S ax and the 
dominant contributions in Eq. (44) come from eigenstates 
with energy Ej = E x + wo in the ” sum and from eigen¬ 
states with energy E 1 = E x — wo in the “+” sum. Therefore 
Eq. (44) is well approximated by 

d p (t) = £atcos(ujot) E± E 1^217 I 2 . (45) 

± 7:f2 7a; =±wo 

As expected the dipole moment oscillates at the same fre¬ 
quency of the electric field. Unlike the spectrum of an ultra- 
short probe, in the monochromatic case S(ui) has at most one 
peak. For x = g (ground state) the ” sum vanishes since 
Cl-yg > 0, and we recover the well know physical interpreta¬ 
tion of equilibrium PA experiments: peaks in S(ui) occur in 
correspondence of the energy of a charge neutral excitation. 
This remains true for x ^ g but the sign of the oscillation 
amplitude can be either positive or negative. Notice that the 
oscillation amplitude is proportional to h x in Eq. (42) and it is 
independent of the delay. 

In the nonstationary case |'T) is a superposition of nonde¬ 
generate eigenstates. For a fixed f3 in Eq. (44) the domi¬ 
nant contributions come from eigenstates with energy E~ = 
Ep + wq in the ” sum and from eigenstates with energy 


E-y = Ep — wo in the “+” sum. Writing VL ai = Cl a p + U/ 3 7 

it is straightforward to show that 

d P {t)=£ 0 t -Y.Y, ±c >^ t 

otfi =h 

x(e» ^ d ai d^ + e ±iUot J2 da'ydryp'j. (46) 

7:^ 7/ 3=±cuo 

As anticipated below Eq. (34) d p (t) is not monochromatic in 
a nonstationary situation, the frequencies of the oscillations 
being wo ± £l a p- Can these extra frequencies be seen in the 
TR-PA spectrum? The answer is affirmative since e(w) is a 
broad function centered in wo and hence for not too 

large Im[e*(wo ± £l a p)d p (u }o ± Q a p)\ is nonvanishing. Fur¬ 
thermore, the induced electric field is sizable and hence the 
second term in the right hand side of Eq. (34) cannot be dis¬ 
carded. 

The dipole moment in Eq. (46) is substantially different 
from the ultrafast probe-induced d p of Eq. (39). In order to 
appreciate the difference we calculate the amplitude of the 
dipole oscillation of frequency wo and compare it with the 
height h in Eq. (43). By restricting the sum over a, (3 to states 
with Cl a p = 0 in Eq. (46), we obtain the harmonics d p ^ 0 ( t ) 
with frequency ±wo 

d P ,UQ (f) — £o t cos(wof) ^ ^ ) T ^ ) d a rydryp. 

ot(3:Q, a p= 0 ± 7 :n 7Q ,=±cuo 

(47) 

The main difference between the amplitude of t cos(wot) in 
Eq. (47) and the peak height in Eq. (43) is that in the former 
we have a constrained sum over a, /3. Consequently no co¬ 
herent oscillations as a function of the delay r are observed in 
the TR-PA spectrum around frequency wo, in agreement with 
recent experimental findings in Ref. 35. 

V. OVERLAPPING PUMP AND PROBE 

In the overlapping regime the difficulty in extracting phys¬ 
ical information from the nonequilibrium response function 
is due to the presence of the pump in the evolution operator. 
Nevertheless, some analytic progress can still be made for the 
relevant case of ultrashort probes. If the pump is active for a 
long enough time before and after the probe then we can ap¬ 
proximate it with an everlasting field. In this Section we study 
the nonequilibrium response function of systems driven out 
of equilibrium by a strong periodic pump field. As we shall 
see in Section VI this analysis will help the interpretation of 
TR-PA spectra. 

Let us consider a periodic Hamiltonian H(t) = // (t + 

Tp) = with wp = 2-jr/Tp. According to 

the Floquet theorem the evolution operator can be expanded 

as 98 

U(t,t') =Ee- i «“ (t - t ' ) |$ a (l))<$ a (t')|. (48) 

a 

In this equation |^ a (f)) = | 'S? a (t+T P )) = einupt 

are the quasi-eigenstates and Q a are the quasi-energies. 


dp(t) - ^ ^ C < yCl ) (i Ct ryd'yfl 


atfl 7 ± 
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They are found by solving the Floquet eigenvalue problem 

Ek^nk]^) = g Q |f'i n) ) with ii F , nk = fc ) + 
nujpS n k . It is easy to show that if { Q a , |'D^)} is a solution 

then { Q' a = Q a + mup, = |T'£ n_m ^)} is a solution 

too. These two solutions, however, are not independent since 
e - l Q°‘ t \'i’ a (t)) = In Eq. (48) the sum over 

a is restricted to independent solutions. For time-independent 
Hamiltonians (fP”t = 0 for all n / 0) the independent solu¬ 
tions reduce to the eigenvalues E a and eigenvectors 'I',,) of 
H = iT (0) . 

We expand the ground state |\P 9 ) = 

Y2 a h a e _I< ^“ to |T , 0 ! (<o)) in quasi-eigenstates and insert 
Eq. (48) into Eq. (36) to derive the following Lehmann-like 
representation of the nonequilibrium response function 

bgbpe lQ °‘~' t e in '' fit ’d ai {t)d 1 p(t')-c.c. 

a.fi 7 

(49) 

In this formula Q a p = Q a — Qp are the quasi-energy dif¬ 
ferences and d a fj(t) = (T' ct (f)|(i|T , ) a(f)) are the time-periodic 
dipole matrix elements in the quasi-eigenstate basis. 

Let us compare the response function during the action of 
the pump, Eq. (49), with the response function after the ac¬ 
tion of the pump, Eq. (38). Unlike the coefficients c a of the 
expansion of |\P) (this is the state of the system after a time 
r from the switch-off time of the pump) the coefficients b,, 
of the expansion of the ground state |T g ) are independent 
of the delay. Bearing this difference in mind we can repeat 
step by step the derivations of Section IV A and IV B with 
E a ->■ Qa and d a p -A d a p{t) = J2 n e mwp *d^. We then 
conclude that the absorption regions occur in correspondence 


of the quasi-energy differences and of their replicas (shifted 
by integer multiples of ujp). 

It would be valuable to relate the quasi-energies to the pe¬ 
riod and intensity of the pump field. In general, however, this 
relation is extremely complicated. In the following we discuss 
the special case of monochromatic pumps, which is also rel¬ 
evant when treating other periodic fields in the rotating wave 
approximation. 99 ~ 102 

A. Monochromatic Pumps 

For a monochromatic pump the time-dependent light- 
matter interaction Hamiltonian has the general form Hp(t) = 
p e -iu>pt _|_ pt e iupt' p| ence fji o) — fj j s ^e Hamiltonian 
of the unperturbed system, P= P, H W = Pt and 
Il lrL> = 0 for all n > 1. Then the Floquet operator reads 



/ 

• pt 

P — up 

P 

0 

0 • 

\ 


P_F = 


• 0 

Pt 

H 

P 

0 • 


. (50) 



• 0 

0 

Pt 

H + up 

P • 




V 






) 



The operator T~Lp acts on the direct sum of infinite Fock 
spaces. The tridiagonal block structure allows for reducing 
the dimensionality of the Floquet eigenvalue problem. With 
a standard embedding technique it is easy to show that the 
quasi-energies Q a and the zero-th harmonic T 0I) ) of the 
quasi-eigenstates are solutions of P e ff(Q) |tlX 0 )) = Q|'T^ 0 )) 
where 


H eS (Q) = P+Pt 


Q-P-wp-Pt 


Q—H—2lj p — P t 


-p 


-p+p- 


Q — H - j- ujp — P 


Q — H — 3uj p — ■ 


Q—H+2ujp—P 


Q-H+3uj p - 


-P t 


P t 


-Pt. 


i51) 


In the large-wp limit the leading contribution is 1 l,p(Q) = 
H + ^P[P,Pl], which can be diagonalized to address the 
high-energy spectral features. 103104 

The Floquet eigenvalue problem simplifies considerably if 
we retain only matrix elements P a p = (T , a |P|T , ) a) withPo,— 
E/s ~ up, and if the subsets of indices a and j3 are disjoint. In 
this case the Fock space can be divided into two subspaces A 
and B with the property that P = YY aF A /3eB 
We write a state |\&} in Fock space as 14/^} + I'I'b} with 
Id/x) = 1^ 4 )’ X = A,B. Similarly, we split 

the total Hamiltonian P = Ha + Hb into the sum of an op¬ 
erator Ha = YlaeA acting on subspace A and 

an operator Hb = Y2peB acting on subspace 

B. It is straightforward to verify that the Floquet eigenvalue 
problem decouples into pairs of equivalent equations involv¬ 


ing two consecutive blocks. Choosing for instance the blocks 
with n = 0 and n = 1 we find 


( Ha P 
\ P^ Hb + up 




(52) 


The 2x2 matrix on the left hand side is known as the 
Rabi operator. From the solutions of Eq. (52) we can con¬ 
struct the full set of quasi-eigenstates according to |4'j(f)) = 

I^ja) + e“ pt |T'cp). 105 Thus the quasi-eigenstates contain 
only a single replica. Notice that in the absence of pump 
fields the solutions are either Q j = Pj, I’f^) = | v kj) and 

|Tfj> = 0 or Q i = Pj + up, |4fj> = 0 and |*g>) = |* £ >. 

The single replica of the quasi-eigenstates reflects into a 
single replica of the time-dependent dipole matrix elements. 
In fact, the given form of the operator P implies that the dipole 
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FIG. 3: (Color online) Absorption energies as a function of the pump 
intensity P £ [0, 0.2] for up = 0 [DC (grey) region] and as a func¬ 
tion of lop £ [0, 0.9] for P = 0.2 [AC region] of a four-level quan¬ 
tum system consisting of the states 1 ,a £ A with energy E\ = 0, 
E a = 21.2, and states b, 2 £ B with energies Ef, = 20.6 E 2 = 23.1 
(energies in eV). The absorption at energy D 12 is independent of 
P and not shown. The color scale indicates the oscillator strength 
of the transition [(vE'i°i|d|'I'!) 1 2) for D Q7 + u P , for 

fla-y — a >p, and the sum of the two for lop = 0, see Eq. (54)] 
normalized to D = (\E , a |d|\£ r l 8) with a = 1, a and /3 = b, 2. In 
the AC region we also display the absorption energy cup with oscil¬ 
lator strength and —cup with oscillator strength 


operator d couples states in subspace A to states in subspace 
B and viceversa. Therefore 

dap(t) = +e- iupt (^lmZ). (53) 

Inserting this result into Eq. (49) we obtain the nonequilibrium 
response function (i > 0) 


zxM) = 

a/?7 

+ d 7/3 (0) - c.c(54) 


(i)t 


For ultra-short probes e(t) = eo S(t) the induced dipole mo¬ 
ment d p (t) = eox(f, 0)- and we can easily deduce the po¬ 
sition of the absorption regions in the TR-PA spectrum from 
Eq. (54). 

We conclude this Section by discussing the paradigmatic 
situation of a pump coupling only two states, say a and b. 


Then P a p = P for a = a and /3 = b, and zero otherwise. The 
quasi-energies are 

E a + E b + up ± ^/{E^~^~~upy~^4P^ 

Q± = -^-, (55) 

and for £ ^ a,b, = E^ for (64 and = E^ + up for 
£ £ B. In the limit of zero pump intensity P —f 0 the quasi¬ 
energies Q + —» E a and —> E b + cup, as it should be. Let 

us analyze with the help of Fig. 3 the various time-dependent 
contributions in Eq. (54). For both a and 7 different from ± 
the square bracket is nonvanishing only provided that a £ A 
[B] and 7 £ B [A]. More precisely, only the first [second] 
term is nonvanishing, and it contributes with an oscillating 
exponential of frequency E a — E 1 — cup + cup = E a — E 1 
[E a + cup — Ery — top = E a — Ej], Thus, the absorp¬ 
tion energy between “pump-invisible” states is preserved (and 
hence not shown in Fig. 3). The situation is more interest¬ 
ing for a £ A [B] a “pump-invisible” state and 7 = ± , 
see Fig. 3 where the “pump-invisible” state is a = 2 £ B. 
Again, only the first [second] term in the square bracket is 
nonvanishing and the corresponding oscillation frequency is 
E a — Q± + up [E a — Q±]. We observe that for P = 0 
the quasi-eigenstate | + (t)) = I'Pq) is in the A subspace 
whereas the quasi-eigenstate |\D_(t)) = e lwpt |4 ; f,} is in the B 
subspace. Therefore the absorption at energy E a — Q + + up 
[E a — Q_ (see the line E 2 — Q- in Fig. 3)] is prohibited 
by the dipole selection rule. The pump field mixes a and 
b, thereby giving rise to the appearance of a new peak for 
every equilibrium-forbidden transition between the “pump- 
invisible” state a £ B [A] and the state b [a]. For cup far from 
the resonant frequency f2 a j, = E a Ei, the (allowed) equilib¬ 
rium transition E a — E^ —> E a — Q_ + cup [E a — E a —> 
E a — Q+], thereby undergoing a shift known as the AC Stark 
shift . 106 - 107 At the resonance frequency cup = the quasi¬ 
energies Q± = E a ± P and the equilibrium peak at E a — Ef, 
[E a — E a ] is replaced by two peaks of equal intensity at en¬ 
ergy E a — Eb T P [ En - E„ + / J ]. This spectral feature 
is known as the Autler-Townes doublet or splitting (since the 
original equilibrium peak appears split in two ). 106108 A sim¬ 
ilar analysis applies for a = ± and 7 £ B [A] a “pump- 
invisible” state. Finally we consider the contributions with 
a = + and 7 = — [a = — and 7 = +] in Eq. (54). In this 
case both terms in the square brackets contribute and the equi¬ 
librium peak at energy Q a b = E a — Eb [D/,„ = Ef, E a ] splits 
into two peaks at energy Q + — Q_ ± cup [Q- — Q+ ± cup], 
see Fig. 3. It is worth noticing that this splitting and the 
Autler-Townes splitting have different origin. In the latter a 
prohibited transition becomes an allowed transition whereas 
in the former a genuinely new transition appear. At the res¬ 
onance frequency the spectrum exhibits two peaks of equal 
intensity at energies Q + — Q_ ± cup = 2 P ± ( E a — Eb) 
[Q~ — Qp±up = —2 P±(E a — Eb)]. Therefore, in this case 
too the equilibrium peak at energy f l a b [fl/ m ] appears split in 
two. Unlike in the Autler-Townes splitting, however, the dis¬ 
tance between the peaks is 4 P instead of 2 P, with the peak 
at energy f+ 2 P [ih m — 2 P] stemming from a shift of the 
equilibrium peak at energy tt a b [Up,], and the peak at energy 
f l a b — 2 P [Af, n + 2 P] stemming from the newly generated 
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transition associated to the equilibrium peak at energy 1 
[(>„/,]. In addition to all the aforementioned absorption fre¬ 
quencies we have the pump frequency. In fact, for a = 7 = =t 
the square bracket in Eq. (54) is the sum of two oscillating 
exponentials with frequency ±wp. Thus, at the resonance 
frequency the spectrum exhibits a three-prong fork structure 
known as the Mollow triplet : 106 109 the side peaks at energy 
f l a b ± 2 P [f 2 h a ± 2 P] and a peak in the middle at energy 

LUp = [ LOp — flfra]. 


VI. MORE ANALYTIC RESULTS AND NUMERICAL 
SIMULATIONS 

The analysis of the nonequilibrium response function x car¬ 
ried out in the previous Section is useful for the physical inter¬ 
pretation of TR-PA spectra. In practice, however, it is numeri¬ 
cally more advantageous to calculate the probe induced dipole 
moment d p directly. In this Section we present some more an¬ 
alytic results for systems consisting of a few levels and single 
out the effects on the TR-PA spectrum of a finite duration of 
the pump. 

Let p be the many-body density matrix in, e.g., the eigen- 
basis of the unperturbed Hamiltonian. The (a, 0) matrix ele¬ 
ment of p is therefore ('I , Q \p\ T^). In the same basis the matrix 
which represents the unperturbed Hamiltonian is diagonal and 
reads H_ = diag( { E a }). We use the convention that an under¬ 
lined quantity 0_ represents the matrix of the operator O in the 
energy eigenbasis. The time-evolution of the density matrix is 
determined by the Liouville equation 

= \H + H P (t) + H p (t), p{t)\ - T,p(t )}, (56) 

where T is a decay-rate matrix accounting for radiative, ion¬ 
ization and other decay-channels. In Eq. (56) the matrices H P 
and H_ p represent the pump and probe interaction Hamiltoni¬ 
ans, and the symbol “[ , ]” (“{ , }”) is a (anti)commutator. 
We set the switch-on time of the pump at t = 0 (hence the 
probe is switched on at time t = r). As we are interested in 
the solution of Eq. (56) to lowest order in the probe field we 
write p = p + p , where p p [t) is the time-dependent den¬ 
sity matrix with H_ p = 0. Then, the probe-induced variation 
P p satisfies (omitting the time argument) 

ij t P p = [H + H P ,p p ] + [H pl p p } - '-{T,p p }, (57) 

which should be solved with boundary condition p (t) = 0. 
For ultra-short probes H_ p (t) = S(t — r)eod and for times 
t > t Eq. (57) simplifies to 

ij t P p = \H + H P ,p p )-^{T,P p }, (58) 

which should be solved with boundary condition p (t) = 
—ieo[d, p p (r)\. Once p p is known the probe-induced dipole 
moment can be calculated by tracing: d p = Tr [p p d\. 



FIG. 4: (Color online) TR-PA spectrum (normalized to the maxi¬ 
mum) for the two level system described in the main text. The pa¬ 
rameters are fl a b = E a — Eb = 7, P = 0.7, 7 a = 7b = 7 = 0.2, 
7 p = 0.02 and u>p = fl a b- Energies in eV and times in fs. 

A. Two-level system 

We consider a two-level system with unpertubed Hamilto¬ 
nian H_ = diag (E a ,Eb), decay-rate matrix T = diag( 7 Q , 7 &), 
and an exponentially decaying pump field which is suddenly 
switched-on at time t = 0 : 

* ( 0 p-iupt \ 

H P (t) = 0{t)Pe e Zt Q J. (59) 

At time t < 0 the state of the system is |T a ) and hence the 
density matrix p (t < 0) = diag(l,0). In Fig. 4 we show 

the TR-PA spectrum S(u>) for 7 p <C 7 , uip = fl a b and for 
an ultra-short probe H_ p (t) = S (t — r)eod, where the dipole 
matrix has off-diagonal elements d a b = db a = ® and zero on 
the diagonal. The spectrum is calculated from Eq. (34) with¬ 
out the inclusion of the quadratic term in d p (thin samples). 
We clearly distinguish a Mollow triplet for small r. As the 
delay increases the side peaks approach the peak at top and 
eventually merge with it. The r-dependent shift of the side 
peaks is a consequence of the finite duration of the pump. In 
the resonant case we have found a simple analytic solution for 
the probe-induced dipole moment (for simplicity we consider 
la=lb = 7 ) 

1 — e~ lpt 

d p (t) = 2 eoS 2 e _7t [cos(wpr) sin(wpi) cos ( 2 P ---) 

7 p 

1 — e~ lpT 

— sin(wpr) cos (2 P ---) cos(wpf)], (60) 

7 p 

for t > t and zero otherwise. As the Fourier transform d p (oj) 
is dominated by the behavior of d p (t) for times r < t < I /7 
we study Eq. (60) in this range. 

We write t = r + s and define the function a(t) = 2P(1 — 
e 1pt )l r ) p. For 7 p <C 7 (this is the situation in Fig. 4) we 
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FIG. 5: (Color online) TR-PA spectrum (normalized to the maxi¬ 
mum) for the two level system described in the main text at a de¬ 
lay r = it/Aup. The parameters are Q a b = E a — Eb = 7, 
7 a = = 7 p = 0.04. Different panels refer to different pump 

intensity P = 0.2, 0.4, 0.6, 0.8 and the solid (dashed) line refers 
to the resonant (off-resonant) pump frequency lop = fl a b + Aw. 
Energies in eV. 


can approximate 

a(t) = a{T) + 2Pe~ lpT s +... (61) 

We then see that for a fixed r the first term in the square 
brackets in Eq. (60) oscillates (as a function of s) at fre¬ 
quencies uip ± 2 Pe~ lpT whereas the second term oscillates 
at frequency up. This explains the shrinkage of the split¬ 
ting between the two side peaks with increasing r shown 
in Fig. 4. Another feature revealed by Eq. (60) is that 
the side-peak intensity is proportional to cos (upr) while the 
central-peak intensity is proportional to sin(wpr) (antiphase). 
The Mollow triplet is therefore best visible only for delays 
r = (2 n + l)7r/(4wp), with n integers. 

Equation (60) is valid for arbitrary 7 and 7 p. Hence, it can 
be used to investigate regimes other than 7 p -C 7 . In the op¬ 
posite regime 7 p > 7 we can approximate a(t) ~ 2P/ r yp, 
and only one peak at frequency up is visible in the spectrum. 
The intermediate regime, 7 p > 7 , is definitely the most inter¬ 
esting as it is characterized by a nontrivial sub-splitting struc¬ 
ture. In Fig. 5 we show the transient spectrum for 7 = 7 p = 
0.04 at delay r = 7 t /(4 up) for different pump strengths P. 
The results are obtained from the solution of Eq. (58). We 
considered a resonant pump frequency, up = tt a b, as well as 
an off-resonant one, up = Q a b + A u. Although an analytic 
formula for d p exists in the non-resonant case too, it is much 
less transparent than Eq. (60) and not worth it to present. After 
a careful study of Eq. (60) we found that the number of peaks 
in the frequency range [— 2 P + up. 2 P + up] grows roughly 
like 0.6 x P/jp and that the peak positions tend to accumu¬ 
late around up (in the limit P/'jp —> 00 the frequency up 
becomes an accumulation point). The same qualitative behav¬ 
ior is observed for an off-resonant pump-frequency, the main 
difference being a shift by Aw of the sub-splitting structure. 
The sub-peaks are probably the most remarkable feature of the 



JL 

0.2 
- 0.2 
J -06 

FIG. 6: (Color online) TR-PA spectrum (normalized to the maxi¬ 
mum) for the two level system described in the main text. The pa¬ 
rameters are f l a b = E a — E b = 7, 7 a = 7 t = 7 p = 0.04, P = 0.4 
and up = flab- Energies in eV and times in fs. 


complicated functional dependence of the TR-PA spectrum on 
the pump and probe fields. The sub-peaks are not related to 
transitions between light-dressed states and they can be ob¬ 
served only for pump pulses of duration comparable with the 
dipole decay time. 

The full transient spectrum in the intermediate regime is 
shown in Fig. 6 . We use the same parameters as in Fig. 5 and 
set the pump intensity P = 0.4. Thus, for r = 7 r/( 4 up) 
the spectrum is identical to the one shown in the top-right 
panel of Fig. 5. The sub-slitting structure evolves similarly 
as the main side-peaks. We observe periodic revivals of the 
sub-peaks whose positions get progressively closer to up and 
whose intensities decrease with r. Another interesting feature 
is that for finite 7 p the broadening of the peaks depends on 
r. It is therefore important to take into account the finite du¬ 
ration of the pump when estimating the excitation life-times 
from the experimental widths. 


B. Three-level system 

We consider a three-level system with unperturbed Hamil¬ 
tonian H_ = diag(£'i, Eb, E a ) and decay-rate matrix T = 
diag( 0 , 7 , 7 ). The system is perturbed by an exponentially 
decaying pump field which is suddenly switched-on at time 
t = 0 and couples levels a and b: 

/ 0 0 0 \ 

H P {t) = 6(t)Pe~ lpt 0 0 e iujpt . (62) 

\ 0 e~ iujpt 0 J 

In the absence of the probe the system is in the ground state 
|\I/i) before the pump is switched on, hence p p (t < 0 ) = 
diag(l, 0, 0). At time r we switch on an ultrashort probe field 
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FIG. 7: (Color online) TR-PA spectrum (normalized to the maxi¬ 
mum) for the three level system described in the main text. The 
parameters are E\ = 0, Eb = 20.6, E a = 21.2, P = 0.7, 7 = 0.2, 
7 p = 0.04 and lop = Q a b (top panel) and cup = Q. a b — 0.5 (bottom 
panel). Energies in eV and time in fs. 

HAt ) = e(t)d where e(t) = eo S(t — r) and, for simplicity, 
we take the dipole matrix of the form 

/ 0 0 £ \ 

d= 0 0 0 (63) 

\® 0 0 / 

(we therefore neglect the matrix element ('T a |^l'£ , &))- For 
uip ~ E a — Eb in the near-infrared region and E a — Ei in 
the extreme ultraviolet a similar model has been considered 
by several authors in studies of TR-PA of He atoms (in this 
case 1 = Is 2 , b = ls2s, a = ls2p represent the first three 
levels of He ), 26 ’ 29,31 ’ 35 ' 70 ’ 72 ’ 110-112 although the time-dependent 
probe and pump fields were different. 

In Fig. 7 we show the TR-PA spectrum for a resonant (top 
panel) and off-resonant (bottom panel) pump frequency cup. 
The spectrum is again calculated from Eq. (34) without the 
inclusion of the quadratic term in d p (thin samples). For neg¬ 
ative r we only see the equilibrium peak at frequency H a i, 
which corresponds to the excitation from the occupied level 1 
to the empty level a (given the choice of the dipole matrix in 


Eq. (63) this is the only possible transition). For small positive 
t we recognize the Autler-Townes splitting discussed in Sec¬ 
tion IV B. The finite duration of the pump causes the collapse 
of the Autler-Townes splitting as r increases; for r —> 00 
we eventually recover the equilibrium PA spectrum. The 
r-dependent shift of the Autler-Townes peaks can be study 
quantitatively in the resonant case. In fact, the differential 
equation for p , see Eq. (58), can be solved analytically for 
uip = il a i, and the corresponding probe-induced dipole mo¬ 
ment reads 

dp(t) = 8(t — T)2eo'£> 2 e~~ ,< ' t ~ T ' > sin — r)) 

1 _ e -lP(t-T) 

x cos (Pe 1PT -) (64) 

7 p 

Interestingly d p {r + s) depends on r only through the expo¬ 
nentially renormalized pump intensity Pe _7pT : the PA spec¬ 
trum at finite r and pump intensity P is the same as the PA 
spectrum at r = 0 and pump intensity Pe~ lpT . There¬ 
fore, the Autler-Townes splitting follows the exponential de¬ 
cay of the pump, in agreement with the numerical simulation 
in Fig. 7. 

In the TR-PA spectrum of Fig. 7 we have 7 p -C 7 . How¬ 
ever, the analytic solution in Eq. (64) is valid for all 7 and 7 p. 
Like in the two-level system a sub-splitting structure emerges 
in the intermediate regime 7 ~ "f P (not shown). A simi¬ 
lar finding was recently found for trigonometric and square 
pump-envelops . 111 

The TR-PA spectra shown so far have been calculated under 
the assumption that the total electric probe field acting on the 
electrons is the same as the external (bare) field. As discussed 
in Section IV B, this approximation makes sense for thin sam¬ 
ples. For samples of thickness much larger than the inverse 
transition energies of interest, the Liouville equation for the 
density matrix, Eq. (56), should be coupled to the equation 
for the total electric field, Eq. (24). 23 ’ 35,40 ’ 70 ’ 75 ’ 90-97 To appre¬ 
ciate the qualitative difference introduced by a self-consistent 
treatment of the probe field we consider again the system of 
Fig. 7 and add to the bare (5-like probe an induced exponen¬ 
tially decaying planewave of frequency 0„i. We enforce (for 
simplicity) monochromaticity on the probe Hamiltonian (or 
equivalently we work in the rotating wave approximation) and 
take a total probe field £ p (t) = (£ x (t — r),£ y (t — t),0), 
where the components 

£ x (t ) = eoS(t) + d(t)\e~ lpt cos Q. a \t 

£ y (t) = 9(t)\e~ lpt sin Pl a it (65) 

decay on the same time scale of the pump field. Choosing the 
dipole components d = ( d x ,d y ,d z ) with 

/ 0 0 $ \ / 0 Oi0\ 

4= 0 0 0 ; d y = I 0 0 0, (66) 

\ $ 0 0 / \ -i® 0 0 / 

the probe Hamiltonian reads 

( 0 0 £ p {t-r) \ 

H(t) = £(t)-d = ®\ 0 0 0 ,(67) 

0 / 
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FIG. 8: (Color online) Linear term in d v of the TR-PA spectrum 
(normalized to the maximum) of Eq. (34) for the three level sys¬ 
tem described in the main text. The probe Hamiltonian is given in 
Eq. (67) with A/eo = 0.01. The rest of the parameters are the same 
as in Fig. 7. The pump frequency is top = f l a b (top panel) and 
ojp = Q, ab — 0.5 (bottom panel). Energies in eV and time in fs. 



FIG. 9: (Color online) Plot of \ui d p (uj)\ 2 (normalized to the max¬ 
imum) for the three level system described in the main text. Same 
parameters as in the top panel of Fig. 8. Energies in eV and time in 
fs. 


According to our calculations, the Autler-Townes splitting and 
the pump intensity decay on the same time-scale, see Eq. (64). 
This implies that a spectral shift can occur only provided that 
the amplitude of the pump field is delay dependent as is, for 
instance, the case in thick samples where the pump is dressed 
by an exponentially decaying dipole moment. 

For thick samples the quadratic term in d p in the TR-PA 
spectrum of Eq. (68) can become relevant. This term is always 
negative and hence it can either suppress a positive peak or 
even turn a positive peak into a negative one. The induced 
dipole moment scales linearly with the sample volume V = 
SL. If we introduce the dipole density per unit volume d p = 
dp/V, then the TR-PA spectrum in Eq. (68) can be rewritten 
as 


with£ p (f) = eo5(t) + #(i)Ae in< * lt_7pt . Our modelling of the 
induced probe field is based on the self-consistent results of 
Ref. 35. We instead assume that the pump field does not need 
to be dressed. 

The formula for the spectrum, Eq. (34), has heed derived for 
linearly polarized probe fields. It is straightforward to show 
that for e = e p n ^e^ the generalization is 


S(u) = -2^Im (jde< n >*(w)4 n) (w)) 

n 


2i r 


uj d p (w) 


(- 68 ) 


where dp 1 ' 1 = • dp"\ In Fig. 8 we display the first term 

(linear in d p ) of the transient PA spectrum of Eq. (68). The 
main difference with the spectrum in Fig. 7 is the appear¬ 
ance of an extra peak at frequency < > i a . It is therefore the 
induced probe field to generate the central peak. Furthermore, 
the height of the central peak increases monotonically (no 
coherent oscillations), in agreement with the results of Sec¬ 
tion IV B. One more remark is about the exponential shrink¬ 
age of the Autler-Townes splitting already observed in Fig. 7. 


S(u) 

V 


-2^Im ( w e<">(w)4">(w)) 

n 


2t tL 
c 


uj dp (uj) 


2 , (69) 


from which we see that the contribution of the last term grows 
linearly with the sample thickness. The quantity |wdp(u;)| 2 is 
shown in Fig. 9 in the resonant case ojp = f } a b- As expected, 
the spectral regions where the linear (top panel of Fig. 8) and 
quadratic (Fig. 9) terms are nonvanishing are the same. Never¬ 
theless, the mathematical structure of the peaks is very differ¬ 
ent, as it should be. In fact, in the absence of damping d p (u>) 
is the sum of Dirac 5-functions and hence its square is the sum 
of Dirac 5-functions squared. 


VII. SUMMARY AND CONCLUSIONS 


We have provided a detailed analysis of the nonequilibrium 
dipole response function x, the fundamental physical quan¬ 
tity to be calculated/simulated for interpreting TR-PA spectra. 
Exact and general properties of x have been elucidated and 
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then related to transient spectral features. In the nonoverlap¬ 
ping regime the height of the absorption peaks are strongly 
affected by the shape of the probe pulse. For ultrashort probes 
the peak heights exhibit quantum beats as a function of the 
delay, a signature of the coherent electron motion in the non¬ 
stationary state created by the pump. As the probe duration 
increases the effects of coherence are progressively washed 
out, and the spectrum is progressively suppressed away from 
the probe frequency. The absorption regions are instead inde¬ 
pendent of the delay and occur in correspondence of the neu¬ 
tral excitation energies (not necessarily involving the ground- 
state) of the equilibrium system. For overlapping pump and 
probe the absorption regions cease to be an intrinsic property 
of the equilibrium system and, more generally, the interpre¬ 
tation of the transient spectrum becomes intricate. Analytic 
results for everlasting periodic pump fields are available and 


at the same time useful for the interpretation of TR-PA spec¬ 
tra. The Lehmann-like representation of \ in terms of light- 
dressed states provide a unifying framework for a variety of 
well known phenomena, e.g., the AC-Stark shift, the Autler- 
Townes splitting, the Mollow triplet, the photon replicas, etc. 

The effects of the finite duration of the pump pulse are dif¬ 
ficult to address in general terms. We have considered the 
two- and three-level systems extensively studied in the liter¬ 
ature and derived an exact analytic expression for the time- 
dependent probe-induced dipole moment. Our solution shows 
that for strong enough pump intensities a rich sub-splitting 
structure emerges, in agreement with the recent theoretical 
findings in Ref. 111. We also find agreement with recent ex¬ 
perimental results on He: for long (induced) probe fields, like 
those occurring in thick samples, the absorption peak at the 
probe frequency does not exhibit coherent oscillations. 35 
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